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Heavy quark production has been for years a closely scrutinized process, 
as it has represented one of the very few instances in which experimental 
measurements and next-to-leading order (NLO) QCD predictions seemed 
to be at variance. Such a disagreement appeared surprising, as one expects 
perturbative QCD to be able to handle heavy quark predictions well, 
due to the (relatively) large scale set by their mass and to their small 
hadronization corrections^. Recently, new theoretical developments and 
better use of higher-quality non-perturbative information have however 
greatly reduced the disagreement in bottom production in pp collisions, to 
the point that it does not appear significant anymore ^. At the same time, 
comparisons for top and charm production appear successful, pointing once 
again to a healthy status for this sector of QCD. 



°* Talk given at IFAE Lecce (Italy), April 2003 

^It should however be noted that, while the large heavy quark mass does indeed make 
accurate calculations possible, at the same time it raises the bar: heavy quark production 
rates become calculable and they are a real prediction of QCD, contrary to the light 
quark case. Recalling that the NLO corrections predict a large increase with respect to 
the leading order rates, it should not come as a surprise if this calculation does not yet 
manage to describe perfectly the measurements. 

^Recent comparisons U to new preliminary bottom data from the Tevatron Run II 
show an even better agreement. I am not reviewing here ep and 77 collisions, where 
discrepancies still seem to exist, but where the accuracy of both data and theoretical 
predictions is also lower. 
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Figure 1: Left: theoretical predictions and uncertainties for the top total cross 
section. Right: comparison with CDF and DO preliminary measurements available 
in April 2003. 



Top - The top quark was discovered at the Tevatron during Run I, and 
it has recently been once again observed in the first Run II data. Its total 
cross section, now measured at the slightly higher centre-of-mass energy 
VS = I960 GeV, can therefore be compared to new predictions. While 
NLO corrections [2] for the total cross section were already available at 
the time of the top discovery, improvements of the last few years include 
the implementation of the resummation of soft-gluon (threshold) effects to 
next-to-leading logarithmic (NLL) accuracy JHJ, and the determination of 
parton distribution function sets (PDFs) with associated errors 1^. Figure H 
displays in the left panel the theoretical predictions and corresponding 
uncertainties fSj for the top total cross section at the Tevatron Run II, i.e. 
pp collisions at = 1960 GeV. For each PDF set the crosses correspond 
to the prediction obtained by setting the renormalization and factorization 
equal to 1/2, 1 or 2 times the central scale. The rectangle around each 
cross represents instead the uncertainty related to the given PDF set. It is 
apparent from the right panel of figure^that the experimental uncertainties 
are nowadays still much larger than the theoretical ones, and that the latter 
are dominated by the PDFs uncertainty, due to our poor knowledge of the 
gluon distribution function at large x. This points to a situation where not 
a better calculation, but rather more accurate experimental inputs will be 
needed in order to perform a more compelling comparison for this observable. 

Bottom - Bottom production in pp collisions at UAl and the Tevatron [B] 
has been for years an example of a possible discrepancy between theory 
and data, by factors of 2 to 3. Many explanations for the discrepancy 
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Figure 2: Left: Fit to e+e moment space data. Right: The new prediction for i?+ 
production at the Tevatron, compared to the experimental data 



were proposed, from conventional ones (like the importance of higher orders 
or small-x contributions) to a more exotic one involving supersymmetric 
particles production 7 . 

A careful analysis of the way the data are extracted and/or theoretical 
predictions are evaluated shows however that some common practices need 
to be revised in the light of today's desired accuracy. Bottom quarks are not 
observed as free quarks, as they hadronize into Hf, hadrons before decaying 
to other bottomless final states. The b — > Hf, transition cannot be described 
by perturbative QCD, and it is usually parameterized by convolving the 
momentum distribution of the quark with a phenomenological function - 
usually extracted from e~^e~ data - which accounts for the degradation of 
its momentum in the hadronization process. The bottom quark not being 
directly observed, there is however no way to determine uniquely its 'frag- 
mentation': when measuring the Hi, cross section we only observe the result 
of both a perturbative (hard gluon emissions) and a non-perturbative (soft 
gluons) degradation of the initial quark momentum. Hence, in the theo- 
retical description the two 'steps' must be properly matched, and no single 
function can be a tool good for all purposes and instances. Data presented at 
the unphysical b quark level, the result of a deconvolution performed on real 
Hf, hadron data, may therefore be tainted by the use of an improper phe- 
nomenological parameterization of the fragmentation effects. Fortunately, 
the CDF Collaboration has also published data for the mesons, real 
observable objects, allowing for a safer comparison between theory and ex- 
periment. 

Ref. |S] implements these considerations by properly matching perturba- 
tive (FONLL, i.e. full massive fixed order calculation to NLO accuracy plus 
resummation of \og{pT/m) terms to NLL accuracy) and non-perturbative 
physics, extracting the relevant (i.e. moments around A*" = 5) experimental 




input from e^e~ data employing the same kind of perturbative calculation 
(and the same parameters) which will then be used to calculate the cross 
section at the Tevatron. The results of this analysis are shown in fig. 121 The 
data and the theoretical predictions are compatible within the uncertainties. 

More theoretical progress in the bottom production sector has taken 
place for the heavy quark fragmentation function. Large-A^ moments are 
not needed to describe hadroproduction, but they are well measured in e~^e~ 
collisions, and have an enhanced sensitivity to non-perturbative contribu- 
tions. All finite order perturbative calculations, and even finite logarithmic 
accuracy resummations, are bound to fail eventually when the endpoint re- 
gion is approached. An 'all-logs' calculation can however be performed in 
the large- /3o limit |2j . Including infinitely many orders in the strong coupling 
this so-called Dressed Gluon Exponentiation (DGE) result is necessarily di- 
vergent, as perturbative QCD series are only asymptotic. The divergence 
manifests itself in terms of poles ( "renormalons" ) in the Borel transform, 
which becomes non-invertible unless a regularization prescription is sup- 
plied. The ambiguity of the regularization procedure is related to the miss- 
ing higher twist terms, whose functional form can therefore now be inferred. 
Hence one is left with a perturbative calculation and with the matching 
power corrections which build up its complementary non-perturbative func- 
tion. The full result can be written as a convolution of the two terms, 
D{x) = D^^{x) D^^{x), neither independent of the other. Being the 
functional form for D^^ suggested by the ambiguity, there is no need to 
resort to a phenomenological model, and only a few parameters have to be 
fixed by using the data. The fragmentation function D{x) resums effects 
at the scale m/N, and it is furthermore universal. It can be convoluted 
with a specific coefficient function in order to describe a given process, e.g. 
e^e~. Figure 01 shows how well the data can be described just by fitting the 
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Figure 4: Left: Theory and data for D'^ production in proton-proton collisions at 
the Tevatron Run II. Right: description of the ALEPH D* fragmentation data at 
LEP with the DGE calculation and the non-perturbative function fitted in 9 to B 
fragmentation, with only the mt — > rric replacement jl3j . 
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to a very limited set of moments. The leading non-perturbative effect, 
shift to the left of the perturbative distribution controlled by eiAQCD, is 
compatible with the expectation that this be a typical hadronic scale, i.e. 
~ 300 MeV. It is also worth mentioning how, fitting also Aqqd in order to 
check for consistency, this does indeed return a value fully compatible with 
other determinations of the strong coupling. 

Charm - A QCD calculation for bottom production should also be able 
to reproduce charm data, given the replacement nib ~^ "T-c and a different 
non-perturbative part, since the quark meson transition is quantitatively 
different. 

To verify this, the same theoretical framework which gives a fairly good 
description of the CDF data at the Tevatron has been checked against 
new charm data from the Tevatron Run II |10j . In this case the de- 
termination of the non-perturbative function was slightly more involved, as 
LEP data were not available for all the states {D^, D^, D* and Ds) mea- 
sured in hadronic collisions. Effective fragmentation functions were there- 
fore constructed from the one for D*, which is accurately measured in e~^e~ , 
modeling the F — > P decays and employing a perturbative QCD model ^21 
to parameterize fragmentation into vector (V) and pseudoscalar (P) mesons 
with one free parameter only. The results for are shown in figure |1] (those 
for and D*^ are qualitatively identical). The FONLL calculation, com- 



plemented by non-perturbative information properly extracted from e~^e~ 
data, is clearly able to describe the CDF data within the uncertainties. 

Last but not least, the resummed calculation of ref. |H| and the form it 
predicts for the non-perturbative power corrections, given in eq. can be 
tested against D* fragmentation data in e^e^ collisions. Figure |1] shows 
in the right panel the ALEPH data and the curve given by the same non- 
perturbative parameters fitted to B meson data and listed in fig. |31 the 
only change having been the replacement of the bottom mass by the charm 
mass in both the perturbative and the non-perturbative components of the 
fragmentation function. The fairly good description of the data, at this stage 
to be taken only as a preliminary indication P!B|, shows that the scaling of 
the non-perturbative effects with the heavy quark mass is indeed correctly 
predicted, and that from the numerical point of view the hadronization 
effects are pretty similar in the B and the D sector, though a more refined 
analysis should certainly be performed. 
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